Asymptotically exact probability distribution for the Sinai model with finite drift 
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We obtain the exact asymptotic result for the disorder-averaged probability distribution function 
for a random walk in a biased Sinai model and show that it is characterized by a creeping behavior 
of the displacement moments with time, (x") ~ t''", where /i < 1 is dimensionless mean drift. 
We employ a method originated in quantum diffusion which is based on the exact mapping of 
the problem to an imaginary-time Schrodinger equation. For nonzero drift such an equation has 
an isolated lowest eigenvalue separated by a gap from quasi-continuous excited states, and the 
eigenstate corresponding to the former governs the long-time asymptotic behavior. 
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The Sinai model [T] is one of the simplest models for 
transport in classical disordered systems where physical 
behavior of various disorder-averaged quantities is inter- 
esting and non-trivial and yet amenable to analytical cal- 
culations. The model describes random walks under the 
influence of a spatially homogeneous thermal white noise 
in the presence of spatial disorder represented by a Gaus- 
sian quenched random-drift field v{x). In one dimen- 
sion (ID) such a disorder strongly suppresses diffusion so 
that the mean random- walk displacement in the absence 
of the mean drift is given by [x^it)) ex In** t P3 [I]- In 
higher dimensions, where the properly generalized Sinai 
model describes various physical applications like hop- 
ping transport in the presence of charged or magnetic 
impurities, dynamics of dislocations in glasses, l//-noise 
etc., diffusion can be either suppressed (sub-diffusion) for 
potential random drifts or enhanced (super-diffusion) for 
solenoidal ones [2H1J . The model also exhibits mesoscopic 
fiuctuations |4] within the ensemble of realizations similar 
to the mesoscopic fluctuations [5] in quantum diffusion in 
the Anderson model. In such a situation, any exact ana- 
lytical result for the Sinai model can be potentially of a 
wide interest and applicability. Naturally, the ID case is 
most promising for finding exact solutions. 

The exact analytic results for the most generic quantity 
in the ID Sinai model, the random walks probability dis- 
tribution function (PDF), has been presented by Kesten 
[B] but only in the limit of zero mean drift while the ex- 
tension of this method to the case of an arbitrary drift 
is not known. Although some exact results for the ID 
Sinai model with an arbitrary drift have been obtained 
- e.g. return probability [2], mean first-passage time [7], 
or persistence [H] - none of the methods used for these 
quantities has been generalized for the PDF. 

In the present publication we close this gap by apply- 
ing the method developed [H] in the context of quantum 
diffusion in the presence of losses (non-Hermiticity) to 
obtain exact asymptotic (long-time) results for the PDF 
of the ID Sinai model in the presence of a finite drift. 

The essence of the method is the following. First 



we obtain a formally exact solution for the Laplace im- 
age of the PDF, Veixjx'), for a given realization of the 
quenched random drift field v{x). We represent this so- 
lution in terms of two linearly independent functions, ^f , 
which obey certain boundary conditions only on the left 
or only on the right boundary of the ID system, respec- 
tively. This allows us to formulate the Langevin-type 
equations for ^f, with x serving as an effective time 
variable (and the Laplace variable e being a parameter). 
Then we use the Furutsu-Novikov |10j technique for av- 
eraging over the quenched disorder, expressing the result 
in terms of the eigenvalues of an auxiliary Fokker-Planck 
equations (FPE). Finally, we map FPE to an equivalent 
imaginary-time Schrodinger equation (SE) and obtain an 
asymptotically exact expression (in the long-time limit, 
i.e. for £ — >■ 0) for the disorder-averaged PDF, P^{x~x'). 

The peculiarity of the solution for nonzero drift is that 
the SE is characterized by a two-scale potential and thus 
has an isolated lowest eigenvalue separated by a gap from 
quasi-continuous excited states. The PDF for the biased 
Sinai model is mostly contributed by the former, with the 
latter giving the main contribution to the PDF for zero 
drift. This explains inevitable difficulties for extending 
standard methods successful for the unbiased Sinai model 
to the model with nonzero drift. 

The ID Sinai model (in dimensionless variables [TT]) is 
described by the following Langevin equation: 

^1 + V{X) + 7]{t) . (1) 



dt 



Here ri{t) is the Gaussian thermal noise with ri(t) 
and ri[t)ri{f) = 
noise averaging 







2S{t ~ t') ( ... denotes the thermal- 
^ is a mean drift velocity due to, e.g., 

a constant drag force and v{x) is a quenched Gaussian 

random-drift field with 



{v{x))=0, {v{x)v{x')) =2S{x-x'). (2) 

After averaging over the thermal noise, the Sinai model 
is described in terms of the Fokker-Planck equation: 



{dt - d,V,)V{x, x'- 1) = 5{x - x')5{t) , 



(3) 



Here V^ = dx ^ H ~ v{x) , and P{x, x'; t) is the PDF in 
a given realization of the quenched random field v{x), 
obeying the boundary conditions VxV{x,x';t)\^^^ — 
at the sample boundaries £± -^ ±00. 

Introducing the Laplace transform of the PDF, 

/■OO 

re{x,x')^ dte-'''V{x,x';t), (4) 

Jo 

we rewrite Eq. (|3| as follows, 

{e-dxVx)r,{x,x')^6ix-x'), 
which can be formally solved by the substitution 

Veix,x') ^ 0{x-x')rt{x,x') + e{x'-x)r^{x,x') , 

V-{x)^--{x') (5) 



V:{x,x') 



ipt{x')dx"-Pe {x') - ipe {x')dx'vt{x') 



Here f'^{x) (with a = ±.) are the eigenfunctions of the 
homogeneous equation which satisfy the boundary con- 
ditions only on the left (<p^) or only on the right (<y5^): 



{e-dxVx)p1{x)={)- Vxvl{x)\ 



0. 



Our main task is to find the full PDF, P{x,x';t) — 
{V {x , x' ; t)) , by performing the ensemble averaging over 
the Gaussian random-drift field (l2|. To this end, we in- 
troduce the shifted logarithmic derivatives of the eigen- 
functions 99^: 

aC(a:)=a,[ln(^-(a;)]-/i-z;(x). (6) 

They obey the following Riccati equations 

dx^^ix) = a [(C)2 ~e]-[^i + v{x)] ^^(x) , (7) 

with one-sided boundary condition for each of them, 
^^{x ~ la-) = 0. Then x plays a role of an effective time 
variable and Eqs. ^ a role of Langevin equations de- 
scribing the 'time' evolution of ^~(a;) and the time 'anti- 
evolution' for ^^{x) (as the boundary condition fixes the 
latter at a maximal value of x). For real and positive e 
the functions ^'^ are positive everywhere as S.~ increases 
from and ^^ decreases to 0, as follows from Eq. ([7|. 

Now we can represent the PDF in Eq. ([5]) in terms of 
^f . First, it follows from the definition ^ that 



Pe{x,x') 



^t{x') + ^r{x') 



(8) 



Secondly, we get rid of the explicit dependence on v{x) 
in Eq. ^ by using Eq. ([t]) where we divide both parts by 
^'^ and integrate over x which allows us to connect the 
integrals of ^"^ and over 1/^°^ resulting in 



K{x,x') 



^t{x')+^^{x') 



exp 



—ere 



dxi 

WW) 



(9) 




FIG. 1. The stochastic variable ^~ {x') depends on realization 
of random drifts for the 'earlier time' (< x) while £,^{x) for 
the 'later' (> a;). Therefore, blocks built of f"*" and ^~ are 
statistically independent. 



Finally, we exponentiate the denominator of the pre- 
exponential factor in Eq. iSh as /„°° e~^'^ +'' My which 
allows us to represent V^ as the product of ^~^ and £,~ 
functions. Since the ( ant i) causality of the Riccati equa- 
tions ([7|, ^~{x') depends, in any realization of the ran- 
dom drift field v{x) only on v{x) with the argument 
X < x' while £,^{x') only on those which the argument 
a; > a;', as illustrated in Fig.[l] This allows us to average 
over the Gaussian quenched disorder, Eq. ^, separately 
for ^* functions - in essence, this is the well-known Fu- 
rutsu - Novikov |IOj technique. Thus we represent the 
averaged P^ functions as follows (with cr = ±): 



P^{x,x')^~-dxJ dyR''iy;x,x')Q~^{y;x'), (10) 




where 
K^{y;x,x') = -dy ( exp 



X 

-y^t{-')Tel^ 

J Ce [Xi) 



'-{y;x') 



-yiti^') 



(11) 



The fully averaged PDF, Pe{x, x'), is expressed in terms 
of P^{x,x'), in the same way as the PDF for a given 
realization, Vi;{x,x'), in terms oiP^{x,x'), Eq. ([5]). 

The causality of the Riccati equations allows us to de- 
rive in a standard way [12] the Fokker-Planck type equa- 
tions for the partial distribution functions Q'^ and i?'^: 



(12a) 



adx'Q" = i^ey-ArdyjQ", 

M'' = yiy + l)dy+y{l + afi); 
adx'R" = (ey - dyM") R" = L" R" 



(12b) 



Equation ( |12a ) is solved with the 'initial' condition for 
x' , Q'^(.y',x' = Ig-) = 1, and the boundary conditions 
for the auxiliary variable y, Q'^{y = 0; x') = 1 and 
Q'^(y— J-oo; x') = 0. The appropriate standard solution 



for x' far from the physical boundaries, i.e. \x' — CI >• 1, 
goes over to the stationary, x'-independent function: 



where 2F1 is the hypergeometric function. To match 
these solutions in the region 1 <C y ^ e~^, we expand 



Q''{y-x')^Q%y) = 



K„ 



V£(i + y) 



Eq. (17a) for J/ ^ 1 up to the leading and the first sub- 



K^(2V^) 



(1 + y)^, (13) 



leading term, which yields 
r(2r)y 



k:+^ 



where Kfj_{z) is the modifi ed B essel function. 

The operator 1/°' in Eq. ( 12b I is Hermitian with respect 



m 



r(K + r)r(i 



r) 



+ (r 



to the scalar product defined with the weight function 
w{y) = il + yr'^+' (14) 



and similarly Eq. (17b) for ey ^ 1, which yields 

iey)--+'-' 



m 



2 sin 27rr 



jey)-^- 

r(l-2r) r(l + 2r) 



(18a) 



(18b) 



and It ha s only positive eigenvalues. ^ The solution to ^he terms y-'^"'' and {ey)'''+'- are the main sub-leading 



Eq. (112bf must be everywhere finite, integrable, mono- ^^^ms in Eqs. (fTsS and (flSbl), respectively, provided that 



tonic with y and satisfy the 'initial' condition R'^{y; x' 
0) = —dyQ°'{y), as follows from the definition (11) 



Thus it can be formally represented as R'' (jj;x,x') = 
Q-(x-x)aL QyQ<^(^y-x') . Far from the physical bound- 
aries, i.e. at \x' — £a\ ^ 1, the solution for Q"{y\x') is 
x'-independent, Eq. (13), so that i?°' becomes transla- 
tionally invariant, R°'{y;x,x') — )■ R'^{y;x — x'). Then 
we choose x' = in Eqs. (10) and ^ thus obtaining 
Pf (x) EE Pf (x, x' = 0) as follows: 



P^ix) = -d, dyQ-^y)e 



-xa L' 



dyQ^v) (15) 



The translational invariance of R'^ ensures a useful re- 
lationship, L'^R1{y,x = 0) = aeQ1{y), which follows 
from substituting d^' — —dx into Eq. ( 12b ) and putting 



after applying dx to the definition (11). Substitut 



ing this into Eq. ( 15 ) allows us to reproduce immediately 



a well-known result [^ for the Laplace transform of the 
return probabihty 11.^ = ^/(O): 



n. 



AyQ--Q° 



K^+i(2Vi)K^_i(2Vi) 
K2(2^) 



For an arbitrary x, it is convenient to represent the inte- 



gral in Eq. ( 15 ) in terms of the eigenvalues of the operator 
L" ^ defined by the eigenvalue equation 



L'^P^(y) = ER"E{y) . 



(16) 



First we solve Eq. ( 16 ), considering E as a, given external 
parameter, in the two limits, y ^ 1 and y ^ s^^. We 
are interested in the long-time asymptotic behavior corre- 
sponding to the Laplace variable e <C 1. Therefore, these 
two limiting solutions should match each other in the 
asymptotically wide region 1 ^ ?; ^ e~^. This matching 
allows us to find the actual eigenvalues E{e). 

The two limiting solutions for an arbitrary E are 

R%{y < e-^) = 2Fi(k + r, k - r; 1; -y) , (17a) 

RUy » 1) = y-"K2,.(2Ve^) , (17b) 



K= 1 



^(Tfi, 



= \lf^'-E, 



(17c) 



Re r < 1/2 - which we will show to hold at < /i < 1 [TB] ■ 
In this case, matching the coefficients in Eqs. (18a) and 



(18b) gives the following eigenvalue equation: 



e'- = 



r2(2r) r(-r - cr^/2) r(l - r + cr^/2) 



r2(-2r)r(r- 0-^/2) r(i 



afi/2) 



(19) 



There is a quasi-continuum set of eigenvalues, corre- 
sponding to the imaginary r, which is given for \r\ <^ 



fj. ^ (when the r.h.s. of Eq. (19) is approximately 1) by 

(20) 



^: 



4 



/ Tinm\' 
\ Ine / 



1,2. 



This expression is valid only for ttt, <C | lne| but the lower 
part of the spectrum is all we need to determine the long- 
distance behavior. Note that for /^ = and \r\ ^ 1, 

1, which leads to the shift 
Substituting the 



the r.h.s. of Eq. ( 19 ) equals 

n + 1/2 in the spectrum (20) 



n 



asymptotic eigenfunctions (17) with these values of E 
into Eq. (15) recovers the known result for /i = [3]: 



P,{x) 



4 ^ (-1)" 
7r£ln2(e)^„^rT| 



exp 



11 {r 



Ine 



However, the main difference between the biased and un- 
biased case is the appearance of an additional bound 
eigenstate of L"*" with E <^ /^^/4, corresponding to a 



real positive value of the parameter r, Eq. ( 17c ), close to 
/i/2. In this case, the l.h.s. of Eq. (19) goes to 0, while 
r(r — /^/2) in the denominator of the r.h.s. is close to 
the pole. Expanding it around the pole and substituting 
r — ptjl elsewhere, we find the lowest eigenvalue: 



i?o(£) = 7r/i[sin7r/ir(/i)] 



(21) 



To illustrate finding the eigenstates, it is useful to map 
the eigenvalue equation ( 16 ) for L"*" into an equivalent 



Schrodinger equation with the imaginary 'time' a;, 

-a,*(x, z) = {-dl + V{z)\ *(x, z) , (22) 

with the help of the substitution 



V -^ viz) = sinh^ - 



2 ' 



R^ 



*(a;, z) 



y{z)l[l + y{z)\'^^'' 






U(z) 



Quasi-continuum 




-J — > 



Bound state 



ln(l/e) 



FIG. 2. The effective potential for the Schrodinger equation. 



The effective potential in Eq. (p2| is given by 



U{z) 



/f2 

4 



I6y{z) 



(l + 2/i)(3 + 2^) 
16[1 + y{z)] 



ey 



and is illustrated in Fig. [2] ft underlines the main differ- 
ence between the biased (/x ^ ) and unbiased (^ = 0) 
cases: in the former case the effective potential is char- 
acterized by two scales which results in the emergence of 
the bound state separated by a gap of order ^^/4 from 
the quasi-continuum. 



As a result, only the lowest eigenvalue of Eq. (21 1 con- 
tributes to the large-scale behavior (at x » 4//z'^) of 
P^[x) while P~ is of no interest: the distribution func- 
tion in this limit is one-sided. Thus in the general expan- 
sion for dyQ^ijj] x), 



dyQ+{y- x) = CO i?+(2/) e-^°^ + ^ ^^ Kiv) e" 



^E+a 



we can keep only the corresponding eigenstate Rg ob- 
tained by substituting E ^ Eq in Eq. pSJ). The coef- 



ficients in the above expansion are found as the appro- 
priate scalar products with the weight functions ( |14[ ). In 
the asymptotic limit of Eq. ( 18 1 we find cq = — /i~^. This 



(23) 



leads, with Eq^e) given for 0<//<l by Eq. ( [2T| ), to 

P+ix)=e-^Eo{e)c-^°^''>'=. 



This is the main result of this paper. As an example it 
gives asymptotically exact expressions for the moments: 



'it)) = 



r(n-H) 



^Mr^(A^) 

sinTT/i 



i'"' 



(24) 



This describes a 'creeping' behavior but the propaga- 
tion is considerably faster than in the unbiased case 



(where, e.g., {x'^{t)) ex ln'*i). The results of Eqs. (23) 



and ( 24 ) have already been known [31 [Tl] for a related 
model described by the master equation on a ID lattice, 
Pn = WnPn-i — Wn+iPm with a broad distribution of 



the hopping probabilities W , which diverges at small W 
as tlj{W) ^ W'^~-^. It has been conjectured, albeit with 
some misgivings [2], that these results might be applica- 
ble to the biased Sinai model. Here we have proved that 



Eq. (23) is asymptotically exact. 



Let us finally note Eq.( 23 ) is valid, with an appropriate 



expression for Eo{e), for any fi. Thus for 1< fj,<2, the 
expansion of Eq in e starts from a linear term, 



E^ie) 



fi-1 



(25) 



This means that in this case there is a constant-velocity 
drift, while the correction term defines the anomalous 
sub-diffusive dispersion but precise form of 6^ is of a rel- 
atively little importance. 
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